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We present a study of the heavy-flavor dynamics in nucleus-nucleus 
collisions. The initial (hard) production of c and b quarks is taken from 
NLO pQCD predictions. The presence of a hot medium (a Quark Gluon 
Plasma described by hydrodynamics) affects the final spectra of open-charm 
(beauty) hadrons and their decay electrons with respect to what found in 
pp collisions. The propagation of c and b quarks in the plasma is based 
on a picture of multiple uncorrelated random collisions, described by a 
relativistic Langevin equation. A microscopic evaluation of the transport 
coefficients is provided within a pQCD approach (with proper resummation 
of medium effects). Results for the final spectra of heavy-flavor hadrons 
and decay-electrons are given, with particular emphasis on Raa and 1)%. 

1. Introduction 

Heavy flavor spectra are an interesting probe of the medium formed in 
heavy- ion collisions: c and b quarks are produced in hard processes dur- 
ing the crossing of the two nuclei and the QCD flavor conservation allows 
to "tag" final-state charm and beauty hadrons, so that the color charge 
and mass of the parent parton propagating in the plasma is known. The 
interaction between the heavy quarks and the medium leads to a modifica- 
tion of the final spectra with respect to the pp case. Different mechanisms 
and approaches were proposed in the literature to account for such an effect: 
e.g. medium-induced gluon-radiation [H[2], scattering mediated by resonant 
states [3], calculations based on the Boltzmann equation [3]. As explained 
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Fig. 1. The c (left panel) and b (right panel) momentum diffusion coefficients 
k t/l(p) resulting from the sum of the soft and hard contributions. The sensitivity 
to the intermediate cutoff \t\* ~ m 2 D is very mild. The curves refer to the temper- 
ature T = 400 MeV, with the coupling g evaluated at the scale a = (3/2)ttT. 



at length in |6] we describe the heavy-quark dynamics through the rela- 
tivistic Langevin equation. This amounts to assume that the heavy-quark 
momentum, crossing the plasma, gets changed due to multiple uncorrelated 
random collisions. For independent studies based on the Langevin equation 
see for instance Refs. 0[8j|9]. 



2. The relativistic Langevin equation 

The relativistic Langevin equation |10] is used to study the propagation 
of c and b quarks in the QGP. The variation of the heavy-quark momentum 
in the time-interval At 

^ = -vD( P )p i +e(t), (2.i) 

is given by the sum of a deterministic friction term and a stochastic noise 
term which is completely determined by its two-point temporal cor- 

relator 

(C(t)¥(t')) = b ij (p)5(t-t'), with b ij {p) = K L ( V )fi? + KT (p)(6 ij -pp>). 

(2.2) 

i (Ap 2 ) 

The latter involves the transport coefficients kt{p) = % At anc ^ k l{p) = 
(Ap 2 ) 

At L , representing the average transverse and longitudinal squared- momentum 
acquired per unit time by the quark due to the collisions in the medium. 

The transport coefficients kj* /l are evaluated according to the procedure 
presented in [3]. We introduce an intermediate cutoff \t\* ^m 2 D (t=(P'—P) 2 ) 
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separating hard and soft scatterings. The contribution of hard collisions 
(|t| > |i|*) is evaluated through a kinetic pQCD calculation of the processes 
Q( p )Qi/i -> Q( p ')<li/i and Q(P)d ~> Q{P')9- On the other hand in soft 
collisions (|i| < the exchanged gluon feels the presence of the plasma. 
A resummation of medium effects is thus required and this is provided by 
the Hard Thermal Loop approximation. The final result is given by the 
sum of the two contributions kt/l(p) = k tJl{p) + K jy£,(p) an( ^ ^ s ex P nc it 
expression can be found in Ref. [5j. In Fig. Q] we display the behavior of 
Kt/l f° r c an d b quarks. The sensitivity to the intermediate cutoff |i|* is 
quite small, supporting the validity of the approach. 

3. Heavy flavor in pp and heavy-ion collisions 

The procedure to generate final-state heavy-flavor spectra can be fac- 
torized into some independent steps. Details can be found in [5]. 

• Initial production of QQ pairs in hard pQCD processes. An initial 
sample of c and b quarks is generated using the POWHEG code [TTlj. 
which implements pQCD at NLO. In the AA case nuclear PDFs are 
employed |12j . Heavy quarks are then distributed in the transverse 
plane according to the nuclear overlap function dN/dx±^TAB(x, y) = 
T J ^(x+b/2,y)TB(x — b/2,y) and given a fc^-broadening depending on 
the crossed thickness of nuclear matter. 

• Langevin evolution in the fireball. In the AA case at the proper-time 
r = \ft 2 — z 2 = To one starts following the stochastic Langevin dynamics 
of the quarks in the plasma. The expansion of the background medium 
is described by ideal/viscous hydrodynamics [T3 | [T4" l [T5]. 

• Hadronization and decays. Heavy quarks are made hadronize (around 
the critical energy density e c in the AA case) using Peterson fragmen- 
tation functions [16J with branching fractions into hadrons taken from 
Refs. |17l I18j. Finally each hadron is forced to decay into electrons 
with PYTHIA [19] , using updated tables of branching ratios [20 1 . 

4. Numerical results 

Our predictions for the cross sections of heavy-flavor electrons in pp 
collisions at RHIC (^/s = 200 GeV) and LHC (^s = 5.5 TeV) are given in 
Fig. [2j PHENIX results are well reproduced. Going from RHIC to LHC one 
observes a sizable increase of the cross section, a hardening of the spectrum 
and the growth of the relative contribution of electrons from b decays. 

For AA collisions effects of the Langevin evolution where studied through 
the nuclear modification factor RAA(PT) = (dN/dpT) AA /{N) co n(dN/dpT) pp 
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Fig. 2. Invariant differential cross section of electrons [(e + + e~)/2], from heavy- 
flavor decay in pp collisions at y/s — 200 GeV (left panel) and y/s = 5.5 TeV (right 
panel). The curves refer to electrons from c or b quarks generated by POWHEG. 
PHENIX data are extremely well reproduced. 



RHIC 9 200 GeV, b=8.44 fm, viscous hydro LHC 9 5.5 TeV, b=8.77 fm (0-90%), viscous hydro 




Fig. 3. The Raa of heavy-flavor hadrons and decay-electrons in minimum-bias AA 
collisions at RHIC (0 — 92% of the total hadronic cross section, left panel) and 
LHC@5.5 TeV (0 — 90%, right panel) for viscous hydrodynamics. 



and the elliptic flow coefficient U2(pt) = (cos(2(/>)) Pt . Here we display results 
obtained using viscous hydrodynamics, for few representative cases among 
the ones explored in [5]. In Fig. [3] we show our findings for the Raa(pt) of 
heavy-flavor hadrons and their decay-electrons at RHIC and LHC@5.5 TeV. 
A general agreement with PHENIX data [22J can be attained. Notice that 
the higher temperatures reached at LHC lead to a larger quenching of the 
spectra; however in the inclusive (e c + e{,) result the effect is partially com- 
pensated by the stronger weight acquired by the b contribution. In Fig. 2] 
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Fig. 4. The nuclear modification factor of heavy-flavor electrons at RHIC (a) and 
LHC@5.5 TeV (b) for viscous hydrodynamics (ro — 1 fm), obtained by integrating 
the electron yields over the indicated momentum ranges, as a function of iV par t . The 
centrality dependence found by the PHENIX experiment [35] is nicely reproduced. 
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Fig. 5. (a) The minimum-bias anisotropy parameter Vi for heavy-flavor electrons 
at RHIC for viscous hydrodynamics. PHENIX data[23 [21] are compared with the 
outcomes of our calculations for electrons from c and b quarks and their combina- 
tion, (b) As in panel (a), but for LHC and two choices of tq. 



the centrality dependence of Raa is displayed. The overall agreement with 
PHENIX data looks satisfactory. 

In Fig. [5] we display our results for V2- Even assuming a very rapid 
thermalization, our findings tend to underestimate the experimental data: 
perturbative cross sections favor small-angle scattering. However in the 
Pt regime covered by the experimental measurements hadronization could 
play an important role; coalescence with quarks of the medium, displaying 
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a sizable elliptic flow, could contribute to enhance the heavy-flavor V2 [3]- 

5. Conclusions 

We studied the dynamics of c and b quarks in the medium produced 
in heavy-ion collisions. In our approach, based on the relativistic Langevin 
equation, the quarks interact with the medium through multiple uncor- 
rected random momentum-exchanges, which would asymptotically drive 
them to thermal equilibrium. The heavy-quark transport coefficients were 
evaluated within a weak-coupling framework. Our procedure leads to re- 
sults in reasonable agreement with the available experimental data, which 
extend up to quite high px- This could even suggest to reconsider the role 
of collisional energy-loss in the suppression of light-hadron spectra. 
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